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Abstract 

The infrared bar-spectrum of a single ammonia molecule encapsulated in nano-cage Ceo fullerene 
molecule is modelled using the site inclusion model successfully applied to analyze spectra of CO2 
isotopologues isolated in rare gas matrix. Calculations show that NH3 can rotate freely on a sphere 
of radius 0.184 A around the site centre of the nano-cage and spin freely about its C3 symmetry 
axis. In the static field inside the cage degenerate V3 and ^4 vibrational modes are blue shifted and 
split. When dynamic coupling with translational motion is considered, the spectral signature of the 
1/2 mode is modified with a higher hindering barrier (2451 cm -1 ), an effective reduced mass (6.569 
g.mol -1 ) and a longer tunneling time (55594 ps) for the fundamental level compared to gas-phase 
values (2047 cm -1 ), (2.563 g.mol -1 ) and (20.85 ps). As a result this mode is red shifted. Moreover, 
simulation shows that the changes in the bar-spectrum of the latter mode can be used to probe the 
temperature of the surrounding media in which fullerene is observed. 

1 Introduction 

Quite recently, " buckyballs" have been observed for the first time in space by Astronomers of University 
of Western Ontario [T] [5] using NASA's Spitzer Space Telescope. They identified spectral signatures of 
fullerenes in a cloud of cosmic dust surrounding a distant star 6500 light years away. The existence of 
fullerenes has thus been confirmed as anticipated after the first observation in laboratory by Kroto et al. 
at Rice University [3J . One may therefore expect observation of molecular species trapped in fullerenes in 
interstellar medium. In this work we investigate on the changes that are induced in the spectral signatures 
of species trapped in fullerene in order to determine how these changes can be used to characterize the 
environment where fullerenes are observed. From inert matrix isolation spectroscopy, one knows that 
infrared (IR) spectra of molecular species trapped in a site are simplified in that, as rotational structures 
are absent in the spectra, a vibrational transition is observed as one peak blue or red shifted or sometimes 
two or more in case of splitting of degenerate levels [H [5] [3] . 

Trapping of atoms or ions and small diatomic molecules in fullerene have been studied both exper- 
imentally and theoretically [3 [51 IE HOI HH H21 H3 IH1 HS1 US HZ1 UHl US EDI However, except 
for CO [T71 [T8] [19] and H2 [20] [21] molecules, only a few of the theoretical work has been done in the 
infrared region. In the case of NH3 trapped in fullerene, calculations have been performed from ab-initio 
methods to study the trapping mechanism [55] [53] [53] but no experimental work has been published. 
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Erkoc and Tiirker [23J have reported that up to six NH3 molecules can be trapped in C^o, on one hand 
and Ganji et al. [23] suggest that only one NH3 molecule can form a stable complex NH3@Cgo on the 
other hand. While Slanina et al. [22) have reported a binding energy of -1830 cm -1 with the ammonia 
molecule oriented towards a pair of parallel pentagons. 

The aim of this paper is to determine how fullerene modifies the spectral signature of ammonia 
molecule, in particular in the region of the vibration- inversion mode (namely umbrella mode). 

A theoretical study of an ammonia molecule trapped in fullerene Cgo is thus performed in order to 
simulate its infrared spectrum in the vibration- inversion frequency region. The choice of NH3 is motivated 
by the fact that it contains a nitrogen atom. CO2 could otherwise have been chosen for the simulation. 

To simulate the infrared spectra of NH3 in fullerene, we use the site inclusion model successfully 
applied to analyze spectra of CO2 isotopologues [B] isolated in rare gas matrix. When the molecule is 
trapped in the nano-cage of fullerene, the interaction with carbon atoms in the short range scale modifies 
the translational and rotational degrees of freedom because the molecule is no longer free to move. As 
a result, the umbrella mode which is characterized by the inversion doubling through tunneling of the 
nitrogen atom across the plane defined by the three H atoms is perturbed as discussed below. 

To calculate the bar-spectra of NH 3 trapped in Cgo nano-cage a renormalization procedure is applied 
on the total Hamiltonian of the system. The dynamical coupling between the molecular degrees of freedom 
and those of Cgo nano-cage is then small enough, at least to a first approximation, to justify an initial 
chaos hypothesis. As a result, the density matrix operator of the system can be factorized into a product 
consisting of two terms, one related to the NH3 renormalized optical states and the other to the Cqo bath 
states spanned by its vibrational modes. The dynamical line-shifts and line-widths of the bar-spectra can 
be disregarded and Cgo nano-cage considered as rigid. 

The interaction potential model used in our calculations is presented in section 2 and the renor- 
malization method applied to separate the different motions from each other described. Results of the 
calculation to determine the frequency shifts of the vibrational modes and the orientational level schemes 
of the molecule are then given. Sections 3 and 4 are devoted to a presentation of the formalism used for the 
construction of the infrared spectra of NH3 in fullerene and to a discussion of the infrared bar-spectrum 
in the frequency region of the vi mode. 

2 The interaction potential energy 
2.1 Potential energy model 

The interaction potential energy Vmc between the trapped ammonia molecule NH3 and the rigid fullerene 
molecule Ceo is modelled as the sum of 12-6 Lennard- Jones (LJ) pairwise atom-atom potentials charac- 
terizing the repulsion-dispersion contributions and the induction part due to the interaction between the 
permanent electric multipole moments of the molecule and their images created at the positions of the 
polarized fullerene carbon atoms. It can be expressed as 

60 4 ( , v 12 / ^6") ,60 

where i and j denote the ith atom of the ammonia molecule and the jth carbon atom of the fullerene 
molecule, respectively; ey and cry are the mixed LJ potential parameters, obtained from the usual 
Lorentz-Bcrthelot combination rules ey = y / eI7ejJ and 2cry = an + cry and ry is the distance vector 
between the ith atom of the molecule and the jth carbon atom. 

In the second term of Eq. {1} is the field generated by the molecular permanent electric multipole 
moments (//, 0,...) of the NH3 molecule at position ro on the jth carbon atom of the fullerene molecule 
at position Yj with polarizability tensor a J . It is expressed as 

Ei ' = vv (^)'" + 5 vvv (kF^) :e + - w 



2 



Note that the multipole moments of NH3 and the polarizability tensor of each carbon atom of Ceo 
are usually defined and given with respect to local frames with their origins at the centre of mass of 
NH3 (G,x,y,z) or each carbon atom (Cj,x.j,yj,Zj) respectively. In the calculation, it is then necessary to 
express these quantities in an absolute frame (0,X,Y,Z) defined with respect to the fixed Ceo molecule. 
We use a transformation rotation matrix following Rose's convention [25) and given in Appendix B. 

Figure 1 indicates how to describe the internal degrees of freedom of the molecule with respect to its 
frame (G,x,y,z) and its external orientational and translational degrees of freedom with respect to the 
absolute frame, and Table 1 gives the various fullerene and ammonia characteristic values used in our 
calculations. 




Figure 1: Geometrical characteristics of a molecule trapped in the fullerene molecule. (0,X,Y,Z) and 
(G,x,y,z) represent the absolute frame and the molecular frame, respectively, q and /3 define the equi- 
librium internal coordinates of the trapped molecule. 

The distance vector in Eq.((T|) can be expressed in terms of the position vectors ro of the molecular 
centre of mass (cm.) and Tj of the jth carbon atom of the fullerene molecule with respect to the 
absolute frame (0,X,Y,Z) and of the position vector rji of the ith atom of the molecule with respect to 
its associated frame (G,x,y,z): 

rij = Vj -r -r)i, (3) 

The IR spectra of the trapped NH3 molecule are usually given in terms of molecular vibrational {Q} 
degrees of freedom. 

Thus, to determine the influence of the surroundings on the molecular internal motions, the position 
vectors r\i of the atoms in the molecule, the molecular dipole moment vector fi, and the quadrupole 



3 



Table 1: Numerical parameters for NH3 molecule and for Cgo molecule used in our calculations. 

NH3 C 6 Q 

q (A) 1.0156 RW (A) 3.556 

[3 (deg.) 68. d(C-C)^ (A) 1.455 

d(C=C)W (A) 1.398 

^ c (D) 1.476 (A 3 ) 1.44 

6 C (DA) -2.930 ajj (A 3 ) 0.41 

B (cm" 1 ) 9.941 



C (cm" 1 ) 


6.309 








C - C w 


N - N^ 


H - H (c ) 


e (cm -1 ) 


19.4 


26.5 


17.2 


a (A) 


3.40 


3.38 


2.53 



(a) From Rcf. [26 . (b) From Ref. [27 . (c) From Ref. [36]. 



moment tensor can all be expressed in a series expansion with respect to the molecular frame in terms 
of the vibrational normal coordinates {Q}. 
Taking only the first order terms gives 

m = vt + ^2*iQv, 

v 

= ® e + ^c v Q v . (4) 

V 

In these expressions 77° is the vector position of the ?th atom in the rigid molecule, and /i e and C 
are its permanent multipole moments, while a£, b", and c v are the first derivatives of rji, fi, and © with 
respect to the normal coordinate Q v which describes the z^th molecular vibrational mode with frequency 
oj v . Note that the expression of [i will also be used in the near and far infrared spectra calculations. 

However, it must be noticed that because of the dual nature of the V2 mode of the ammonia molecule, 
that is, the high frequency vibration-low frequency inversion, this mode will be treated in a particular 
way. 

Using the Born-Oppenheimer approximation (adiabatic approximation) to separate the high frequency 
vibrational modes {Q} of the molecule from its low frequency external modes, that is, the orientation 
f2 = ((/j, 0, \) arL d cm. translation ro motions, the interaction potential energy Vmc can be written as 

V MC = V M (r , n) + W M ({Q}) + AU M (r , n, {Q}), (5) 

where Vm(i"o, fi) represents the low frequency motions dependent part of the potential energy for the non 
vibrating molecule, and Wm({<3}) characterizes the vibrational dependent part for the molecule at its 
equilibrium position and orientation, and is generally a small perturbation which involves only vibrational 
level shifts and/or splittings. It can be developed in a Taylor series expansion up to second order with 
respect to the vibrational normal coordinates {Q} as 

™ = £^ + iE^«^ + --- (6 ) 

v vv' 

Finally the third term in Eq.Q characterizes the coupling between the external and the internal 
modes which can lead to the vibrational energy relaxation. 
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2.2 Potential energy surfaces 

2.2.1 Orientation-translation motions 

The first step of the numerical calculations consists in determining the equilibrium configuration of the 
rigid ammonia molecule into the Ceo cage. The potential energy Vm (see Eq.© is minimized with respect 
to both the cm. displacement vector ro and the molecular orientional coordinates (p, 9, and x- 

We find that the equilibrium configuration has an energy minimum value of - 1460 cm -1 for the z 
molecular axis (threefold symetry axis C3) oriented along the direction connected to each two symmetrical 
carbon atoms. The cm. is then displaced from the site centre by about 0.184 A in the same direction. 
An energy maximum of - 1445 cm -1 is obtained when the molecule is oriented and its cm. displaced by 
0.188 A along the pentagon centre directions. There is also an intermediate energy value of - 1455 cm -1 
that is calculated for the molecule oriented and its cm. displaced by 0.180 A along the hexagon centre 
directions. 

Thus, the calculations show that i) in all the molecular orientational configurations, the cm. dis- 
placements ro have always opposite directions to the molecular z axis, and ii) the energy changes are 
negligibly small < 3 cm -1 when ro changes by about 0.01 A along these directions. So, in the following 
we will take the value rg = 0.184 A as a mean value for every orientation. This means that the trajectory 
of the molecular cm. can be approximated by the surface of a sphere of radius r^. 

Note, moreover, that the induction and repulsion-dispersion contributions represent, respectively, 35% 
(- 498 cm -1 ) and 65% (- 963 cm -1 ) of the total energy. 

As a result, the potential energy surface Vm(»"o(^)) is nearly flat since the barrier height is always < 
16 cm -1 and the cm. displacement vectors ro from the cage centre remain always nearly anticollinear 
to the molecular z axis i.e. z- ro ~ - ro. This is not surprising because of the repulsive character of 
the hydrogen atoms of the molecule and reflects a very strong coupling between its orientation and the 
location of its centre of mass. 

Thus, on one hand, the axis of the molecule undergoes only oscillatory motions (ip, 9) of small am- 
plitude around the direction of each displacement vector r{j. The corresponding potential energy surface 
which is given in Figure 2 is a nearly harmonic two dimension librational oscillator. The associated 
frequencies are uj v ~ 155 cm -1 and tug ~ 151 cm . 

On the other hand, for a given molecular orientation, the radial potential function Vm^o) of the 
translation motion of the molecular cm. around the equilibrium position rg = 0.184 A is plotted in 
Figure 3. The curve appears as that of a nearly harmonic one dimensional oscillator with a force constant 
value k = 2.4xl0 4 cm -1 . A corresponding to a frequency w ra( j ~ 218 cm -1 . 

On the basis of these features, we can conclude that ammonia should exhibit a nearly free rotation 
motion simultaneously around the centre of the cage by describing a spherical curve with radius — 
0.184 A, and around its centre of mass. Note however that the spinning motion x remains around the 
threefold symmetry axis of the molecule, only. 

Finally, it is clear that this rotation-cm. translation motion is slow with respect to the librational 
and the radial motions and can be separated from the latter using the adiabatic approximation. 

2.2.2 Inversion-translation motions 

It is also clear, from the results above, that when the inversion motion of the molecule proceeds, its centre 
of mass moves simultaneously from one position (rjjj, SI) to the opposite one (-rjj, SI) with respect to the 
cage centre where the molecule lies in its planar configuration. 

2.2.3 Few remarks 

It must be noticed that the interaction potential energy experienced by the ammonia molecule trapped 
in the fullerene molecule strongly depends on the geometrical characteristics such as the radius of the 
Ceo cage and the equilibrium internal coordinates (q,/3 see Fig. 1) of the ammonia molecule, and also on 
the dispersion-repulsion potential parameters especially the Lennard- Jones a parameters. 

For instance, in one hand, the minimum depth is enhanced from - 1460 to - 750 cm -1 when the 
fullerene radius is decreased by 1%, from 3.556 to 3.519 A, and it diminishes to - 2000 and - 2411 cm -1 
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Figure 2: Potential energy surface versus (f,0) angular motions around the molecular equilibrium con- 
figuration (r§,f2). 

when the radius value is increased by 1% and 2%, to 3.592 and 3.627 A, respectively. The induction 
energy contribution undergoing only weakly variations. 

In the other hand, changes of the a parameters of + 1% and - 1% involve, respectively, minimum 
depths of - 952 and - 1976 cm" 1 . 

However, in spite of these changes in the well depth values, the potential energy surface Vm(^q,^) 
remains nearly flat since its barrier height does not exceed 24 cm -1 . Furthermore, the potential energy 
surfaces associated with the various motions remain nearly unchanged. 

Thus, from the previous analysis of the potential energy experienced by the ammonia molecule into the 
fullcrcnc molecule, we will assume that the cm. position vector ro remains permanently anticollinear to 
the z molecular axis and separate it into a static position vector rg((^, 9, Q2) which parametrically depends 
on the orientation and inversion degrees of freedom, and a dynamical coordinate uo characterizing the 
cm. vibration around (tp, 6, Q2)- 

Finally, it is clear that inversion-cm. translation motions (tq,^), on one hand, and orientation- 
cm. translation motions (rg,!!), on the other hand, are strongly coupled; their kinetic operator must be 
classically treated, as explained in appendix A, before their quantum mechanical treatment as described 
below. The interaction potential energy Vmc can be rewritten as 

V MC = V M (r e , n) + Vm(^o) + V M (r , Q 2 ) + W M ({Q}) + AV M (r Q ,£l, {Q}), (7) 

3 Quantum-classical model 

As we have mentioned above, the purpose of this work is to determine the infrared spectra of the NH3 
molecule trapped in a rigid Cgo molecule. The Hamiltonian of the optically active molecule considered as 
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Figure 3: Potential energy function versus molecular cm. displacement around its equilibrium position 
vector Tq. 

quantum mechanical system can be expressed from the potential energy surfaces calculated above and the 
quantum mechanical model for the orientational and vibrational kinetic energy operators. The latter are 
obtained from the correspondence principle applied to the classical formulation described in Appendix 
A. Disregarding the dynamical couplings between the orientational and vibrational degrees of freedom 
the Hamiltonian can thus be written as 

rreff p/eff _i_ freff . /o\ 

-"a orient ~t~ vib' \°) 

where H°f ient and H^f h are, respectively, the renormalized Hamiltonians associated with the orientation 
and the vibration modes which account for the couplings (rg,f2) an d i r o,Q2)- 

3.1 The orientation modes 

Following expressions derived in Appendix A, the quantum mechanical renormalized Hamiltonian H°f icnt 
for the non- vibrating molecule can be written as 

Orient =2^ + V M (r5,n); (9) 
where T°^ t is the effective rotational kinectic operator 



d 2 d Id 2 ( 2 „ C \ d 2 cote d 2 1 

1- cot ft 1 = 1- cat 2 ft A 9. S- 



T° " = -B cn i — + cot 6— + ^-tt-J + cot^ ft + —s tt-7 - 2 . - 5 5 } . (10) 
rot [dft 2 30 sin 2 ft dip 2 \ B eS J d\ 2 sin 6 depdx) 

In this equation B cS = h 2 /2/g ff is the effective rotational constant connected to the rotation-cm. 
translation motion of the C3 symmetry axis for the trapped molecule and C its rotational constant of 
the spinning motion around this axis for the gas-phase molecule C — 6.309 cm -1 . The value of B is 
equal to 7.419 cm -1 which is less than the gas-phase value of 9.941 cm -1 . 



7 



Furthermore, since the barrier height associated with the orientation-cm. translation motion is 
always < 16 cm -1 , the potential VM(rj),n) can be neglected. The renormalized Hamiltonian becomes 



The eigensolutions of this Hamiltonian are then Ejmk = B eB J(J+l)+(C-B eS )K 2 for the eigenvalues 
(energy levels) and \ JMK) for the eigenvectors where J, M, and K are the usual quantum numbers 
connected to the free rotation motion of the ammonia molecule with the conditions ~J< M < + J and 
— J < K < +J. Note that these levels have (2J+l)-fold degeneracy on the quantum number M and 
twofold degeneracy (±) on the K one, except when K — 29 j. The corresponding line structure of the 
far infrared (FIR) bar-spectrum is then identical to that of NH3 gas-phase one. 

3.2 The vi vibration-inversion mode 
3.2.1 Gas-phase NH 3 

In the vibration-inversion mode of the NH3 molecule, the nitrogen atom can pass from one side of the 
plane defined by the three hydrogen atoms to the opposite one. This mode is generally described, to 
a good approximation, as one-dimension tunneling motion of one particle of mass equal to the reduced 
mass /i2 of NH3 (2.563 g.mol -1 ) moving in a symmetrical double-well potential function V V2 (Q2) (left (L) 
and right (R) wells) as shown in Figure 4a, with a high but finite hindering barrier of 2047 cm -1 |30] , 
such as 



where k 2 = 58306 cm" 1 . A 2 , A 2 = 12469 cm" 1 and a 2 = 4.8224 A 2 . 

The solution of such a problem leads to vibrational energy levels split into doublets, (+) and (— ), 
due to inversion. The vibration- inversion wavefunctions are then described by |i4 ) ( a ~ +> — )j instead 
of \v 2 ), and are expressed as linear combinations (symmetric and antisymmetric) of the vibrational 
wavefunctions \v^) and \v^) associated with the left and right wells, respectively. Moreover, the only 
permitted transitions are (+) < — 5- (— ) as shown in Fig. 4a. 

Note that the potential parameters given above were adjusted to give the fundamental frequency and 
the level splittings corresponding to the experimental values, that is, 950 cm -1 for the frequency and 0.8 
cm -1 for v 2 = and 35.8 cm -1 for v 2 = 1 level splittings. 

3.2.2 Trapped NH 3 

When the NH 3 molecule is trapped in Cgo cage, the gas-phase double- well potential function is modified. 
In the calculations, the inversion-cm. translation potential energy Vyi{r , Q 2 ) (see Eq. 0) is numerically 
added to the potential function V V2 (Q 2 ) of Eq. (fTTj). The new double- well potential function is shown in 
Fig. 4b as a function of the inversion angular coordinate. Its symmetry is preserved but the hindering 
barrier is increased by 404 cm -1 . Moreover, from expressions given in Appendix A, the effective reduced 
mass is equal to fi 2 s — 6.569 g.mol" 1 leading to an increase of 156% with respect to the gas-phase value 
of 2.563 g.mol -1 . 

The corresponding Schrodinger equation was solved numerically using a discrete variable representa- 
tion method |37j . The computed level scheme is given in Fig. 4b. 

In gas-phase, the V2 mode is measured at 931.7 cm" 1 and 968.3 cm" 1 . In Ceo, it is shifted by 47.0 
cm" 1 (as given below) and because the double- well inversion potential is modified with an increase in the 
hindering barrier height and the effective reduced mass as discussed above, this mode is quasi-degenerate 
and is calculated at 751.6 cm -1 , that is, (704.6+47.0) cm -1 (see Fig. 4b). This result is to be related 
to an increase in the tunneling time for both the fundamental and the first excited levels in Ceo, since 
the values are equal to 20.85 ps and 0.47 ps in gas-phase compared to 55594 ps and 333.56 ps in Ceo, 
respectively. We can infer from these results that in Ceo, localization of the N atom occurs in one of the 
wells of the potential function and that the V2 mode is quasi-degenerate. 
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Figure 4: Vibration-inversion double-well potential energy of NH3 as a function of the angle between 
the N— H bond and the plane of the H atoms. The obtained energy levels and the allowed fundamental 
transitions are also shown (a) for the gas-phase and (b) for the trapped one in the fullerene nano-cage. 



3.3 The vibrational frequency shifts 

The second term of Eq. ([8]) is the renormalized vibrational Hamiltonian which accounts for the vibrational 
dependent part Wm({<5}) and inversion-cm. translation part Vm{to,Q2) of the interaction potential 
energy 

Kl = #vib + W M ({Q}) + V M {r ,Q 2 ); (12) 
where -ff v ib characterizes the gas-phase molecular vibrational Hamiltonian written as 

#vib = £|f + ^E KQl + AKTb h ({Q})- (13) 

In this expression /i„ and k v are, respectively, the reduced mass and the harmonic force constant 
connected to the ^th vibrational mode, with associated normal coordinate Q v and conjugate momentum 
P v , and AT^Jjk ({Q}) represents the anharmonic part of the internal potential function. The eigenelements 
of this Hamiltonian -ff v ib have previously been obtained from the ab initio calculations of Martin et al. 
[28]. 

Let E v „ and |. ..«„...) be the eigenvalues and eigenvectors of the vth level associated with the vth 
vibrational mode. 

Thus, a first order perturbation treatment allows us to determine the frequency shift Aw v associated 
with the vibrational fundamental transition | ...0„...) — >| ...1„...) for the trapped molecule from the 
equation 

Aw„ = ft" 1 [{...l u ... I W M {{Q}) I ...1,,...) - (...0,... I W M {{Q}) I ...0,...)] ; (14) 

in which all other modes remain in their fundamental states and h = hjlu. The computed values are given 
in Table 3. We note that each of the vibrational modes are blue shifted with splitting of the degenerate 
levels of v-j, mode (strong) and mode (weak) modelled from static effect. Whereas for degenerate V2 
mode the blue shift of 47.0 cm -1 from static effect, is modified to red shift when dynamic effect on ID 
tunneling motion of the N atom is considered as discussed above. 
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Table 2: Vibrational frequency shifts (cm 1 ) of NH 3 molecule trapped in the Ceo cage, together with the 



gas-phase frequencies. 



Vibrational modes 


Frequency shifts (cm 1 ) 


Gas-phase frequencies (cm L ) 


V\ 


133.0 


3337 


Vl 


47.0 


950 


VZa 


171.0 


3448 


V-ib 


69.6 




VAa 


47.1 


1627 


Vib 


34.0 





4 Infrared absorption spectra 
4.1 General 

The infrared absorption coefficient for the optically active ammonia molecule trapped in a Ceo fullerene 
molecule at temperature T, is defined as the real part of the spectral density, i.e., the Fourier transform 
of the time-dependent autocorrelation function <p(t) 



o 

4>{t) = Tv[p(0)n A (0)n A (t)}, (15) 

where uj is a running frequency variable and c the speed of light, p(0) characterizes the initial canonical 
density matrix operator of the system, and /j,a is the molecular dipole moment operator defined in the 
absolute frame (see Appendix B). The trace (Tr) operation means an average over the initial conditions 
at time t = and over all the possible evolutions of the system between times and t. 

However, since this work is devoted to constructing the infrared bar-spectra of the trapped ammonia 
molecule, we assume the initial chaos hypothesis for the density matrix to be valid. This allows us to 
write the optical wave functions | v v JMK) as products of the renormalized vibrational and orientational 
ones | v v ) <g> \ JMK) Vv . It must be noticed that the orientational states parametrically depend on the 
vibrational states through the moments of inertia. This dependence will be ignored below. 

4.2 Near infrared bar-spectrum 

Within this framework, the infrared bar-spectrum connected to the vth. fundamental vibrational mode 
transition | 0„) — >\ l v ) (all other modes being in their fundamental states) is written in the Lorentzian 
form as 

S~K 2 M ,,„ , ^ ,„ ll2 e -PE jMK _ e -PE 1JlM , K , dp A " 2 



iv{u>) = ^- u \(o v \Q v \i v )\ 2 E 

JMK J'M'K 



JMK 



J'M'K 



dQ v 

X5(C0- UJQJMK^l.J'M'K' - Aw„) . (16) 



In this equation the (...) brackets refer to vibrational transition elements of the normal coordinate Q v 
and to the orientational transition elements of the first derivative of the molecular dipole moment with 
respect to this coordinate. The numerical values connected to all the NH 3 modes are given in Appendix 
B. Moreover, Z defines the orientation canonical partition function associated with the fundamental 
vibrational level at temperature T, (3 = (fc B T) _1 where k B is the Boltzmann's constant, 5 is the Dirac's 
function, and ljojmk-^u'M'K' = hr 1 (Eij/m'K 1 — Eqjmk) is the transition frequency. Note that we use 
below the usual spectroscopic nomenclature to designate the line transitions connected to the rotational 
selection rules, that is, Q(Jk) for AJ = 0, R(Jk) for AJ = +1, and P(Jk) for AJ = —1. 
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4.3 v-i mode results and discussion 



In Figure 5, are given the near infrared bar-spectra of NH3 trapped in fullcrcnc in the V2 mode absorption 
region for three different temperatures, 10 K, 30 K and 100 K. The spectra are calculated by taking into 
account the frequency shift and allowed rovibrational transition elements following selection rules. 

For the latter, as the first derivative of the molecular dipole moment with respect to the Q2 coordinate 
does not depend on the rotational spinning angle Xi = for the K quantum number. 

The infrared bar-spectrum consists of a Q line located at the pure vibration-inversion frequency 751.6 
cm -1 , and a set of rotation-vibration-inversion lines P and R 14.8 cm -1 apart, corresponding to 2B cS , 
on both sides of the Q branch. 




650 
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750 800 
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Figure 5: Near infrared bar-spectrum in the 1/2 vibration- inversion frequency region for NH 3 trapped in 
the fullerene nano-cage at temperatures 10, 30 and 100 K. 

At T = 10 K, only the lowest few rotational states are populated. As a result, the bar-spectrum shown 
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in Figure 5a exhibits only three lines: i) a Q line at 751.6 cm -1 which is in fact a superposition of the 
three transitions Q(li), Q(2i) and Q(22), ii) a R(Oo) line at 766.4 cm -1 and Hi) a R line, superposition 
of R(lo) and R(li) lines at 781.3 cm -1 . There is only one P(lo) line located at 736.8 cm -1 . 

As the temperature increases, more and more rotational states are populated and as result, more lines 
appear in the spectra with a different distribution in terms of intensities. At 30 K, the Q line is the 
most intense consisting of a superposition of Q(Jk) lines for J = 1 to 4 and 1 < K < J. Moreover, the 
intensity of lines R(0o) and R(1at) decreases leading to the increase of the R(2x) line which is located at 
796.1 cm -1 . More lines appear in the R and P branches (R(3k) and P(2k))- 

At 100 K, the Q branch consists of a superposition of 6 lines for J varying from 1 to 6. The R and 
P branches are broader with lines of weaker intensities than those at 10 and 30 K in particular for the 
maxima. 

The modification of the spectra near vi with temperature shows that the spectral signature of the 
trapped molecule can be used to sense the environmental temperature of fullerene Cgo- As the coupling 
with the vibrational modes of Cgo is weak, the signal to noise ratio may be high enough for all the lines 
to be observed. In our opinion this result can be extrapolated to other types of simple molecules trapped 
in C 60 . 

Finally, the dynamical coupling of the vibration-rotation states of the NH 3 molecule with the vibra- 
tional modes of fullerene C 60 involves well separated states for the latter. As a result, this coupling is 
weak and the energy from local translational modes cannot redistribute itself in the thermal bath of Ceo ■ 
The linewidths are then negligible. 

In conclusion, this work shows that the simulation of the spectra of NH 3 can be used to probe the 
temperature of the media in which fullerene is observed. This result can probably be extended to other 
molecules trapped in fullerene, in particular when the coupling with the vibration modes of fullerene is 
weak, because signal to noise ratio should then be high enough for the rotation lines to be observed. 
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A Kinetic Lagrangian expression 

The kinetic Lagrangian associated with the ammonia inversion, orientation, and cm. translation modes 
considered as classical motions is 

Lk = \niQ? + \lB(0 2 + p 2 sin 2 9) + i/ c (^cos0 + x) 2 + ^Mr 2 , (17) 

where /i 2 is the reduced mass of the molecule associated with the Vi mode, Ib and Iq are the molecular 
principal moments of inertia for the rotation motion of its C3 symmetry axis and the spinning motion 
around this axis, respectively, and M is the molecular mass. Moreover, as it has been mentioned above, 
the position vector of the centre of mass of the molecule can be written as ro = ^q(p, 9, Q2) + uq. The 
velocity vector ro is then 

By substituting Eq. (|18p into Eq. (|17p . this latter equation becomes 

1 2 1 1 1 dr e 

Lk = ^fQ 2 + ^If{9 2 + p 2 sin 2 9) + -I G (pcos9 + X f + ^Mu\ + M-^-Q 2 u , (19) 

where /I2 and Ig are effective reduced mass and moment of inertia 
lif = V2+M(J^£\ , and if = I B + Mrf. 

Note that the last term in Eq. (fT9j) accounts for the simultaneous dynamical evolution of the trans- 
lation motion of the centre of mass of the molecule and its inversion motion. 



B Rotational matrix transformation and molecular dipole mo- 
ment 

The unitary matrix M characterizing the transformation from the absolute frame (0,X,Y,Z) into the 
molecular one (G,x,y,z) through the Euler angles ip, 9 and x is given as [25] 

M(p,9, x ) = 

cos p cos 9 cos x — sin <p sin x sin ip cos 9 cos x + cos ip sin x — sin 9 cos x 

— cos (p cos 9 sin x — sin ip cos x ~ sin ip cos 9 sin x + cos (p cos x sin 9 sin \ | . (20) 

cos (p sin 9 sin tp sin 9 cos 9 

Thus, the molecular dipole moment ha and its first derivative vectors §rp, with respect to the normal 
coordinate Q v of the vih vibrational mode can be written in the absolute frame as 

HA = M- 1 ^,^, 
|g = M-\p,9,xW, (21) 
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where M 1 is the inverse matrix of M, and h v — J^- the first derivative vectors of the molecular dipole 
moment expressed in the molecular frame. Their non vanishing component values (in DA ) are [38] 

r bl = fe = i- 404 b2 = % = - 4 - 409 ] 

b3a = lt: = - 4 - 012 b36 = fe = 2 - 465 • 

Moreover, the values of the elements of the normal coordinates connected to the fundamental vibra- 
tional transitions are (in A) 

|(0i|Qi|li)| = 0.042 ] 
|(0 2 |Q 2 |1 2 )| = 0.069 I 
K0 3 |Q 3 |1 3 )| = 0.029 f- 
|(0 4 |Q 4 | 1 4 >| = 0.060 J 
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